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The quantum tunnelling and nucleation theory of vortices in helium II is reviewed. 
\ Arguments are given that the only reliable method to calculate tunnelling prob- 

■ abilities in this highly correlated, strongly interacting many-body system is the 
. semiclassical, large scale approach for evaluation of the tunnelling exponent, which 
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■ represented in some detail and its relevance for the interpretation of experimental 
! data is discussed. 
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1 Preface 



The nucleation theory of quantized vortices in the Bose superfluid helium II 
has been an elusive subject ever since the existence of quantized vortices was 
conjectured by Lars Onsager in 1949. The difficulty to grasp their coming 
into existence in a definite manner from first principles has one fundamental 
reason: There is no microscopic theory of this superfiuid. We do not know how 
to describe the motion of a vortex on scales of the order of the coherence length, 
where this motion is governed by the full quantum many-body structure of the 
superfiuid and the interaction of the vortex with the microscopic excitations. 
Consequently, we can not follow the vortex on any stage of the evolution from a 
virtual vortex to a mature one describable by semiclassical means. At present. 
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even the many-body structure and energy of a simple static rectilinear vortex 
can not be uniquely determined, let alone that of a circular vortex which is 
still beyond reach even for Monte Carlo calculations [1]. 

This lack of a microscopic idea of vortex motion makes it necessary to resort 
to a semiclassical, hydrodynamic (large scale) theory. In such a treatment, 
the existence of the vortex as a semiclassical object has to be assumed ab 
initio. No details of the underlying microscopic dynamics, i.e. of the actual 
nucleation event, are to be described in such a theory, but only the laws 
which rule vortex motion on curvature scales well beyond the core size ^ 
(Fig. 1). This is admissible because in helium II this length scale is of the 
order of the interatomic spacing and the relevant scales for potential barriers 
to be surmounted by thermal activation or crossed by quantum mechanical 
tunnelling usually have spatial extent and length scales far beyond ^. Other 
advantages of ^le include that it can be prepared to have very high purity 
such that dissipation has negligible influence. The order parameter has the 
simplest conceivable structure resulting from a spontaneously broken U(l) 
symmetry. In contrast to the extraordinarily complicated behaviour on (^-scales 
this superfluid can be studied under quite basic and simple assumptions if one 
is remaining completely in the semiclassical, large scale domain. 

I. Quantum Hydrodynamics 
dominated by Magnus force 




Fig. 1. Length scales of vortex motion and interaction with elementary excitations in 
the dense superfluid hcHum II. It is only in the regions I and II where a semiclassical 
picture of vortex motion is applicable. The vortex in the central region of extension 
O(^) is termed virtual because it is no well-defined topological object in this region. 



In this paper, we will discuss and investigate more closely the hydrodynamic 
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theory of vortex quantum tunnelling. In the section which follows we contrast 
it with different other ways to calculate tunnelling rates. Then we explain the 
way in which the classical theory of the tunnelling exponent in three dimen- 
sions is a geometric theory in section 3. How to obtain an explicit solution 
of the boundary problem for a point vortex moving around an ellipse in two 
dimensions is demonstrated in section 4. It is shown that at boundaries with 
strongly varying curvature radii strict bounds on semiclassicality of vortex mo- 
tion are required, which lead in turn to restrictions on the tunnelling energy 
of the vortex. We estimate the prefactor of the tunnelling rate in section 5 and 
the relation to available experimental data is estabhshed. We conclude with 
some general remarks. The discussion will be restricted to zero temperature 
as we will be interested in pure quantum mechanical tunnelling without any 
assistance of thermal fluctuations. 



2 Diflferent approaches 

We present in this section a critical overview of some of the methods developed 
over the last decades to describe vortex quantum nucleation and tunnelling. 

The description of the tunnelling phenomenon by means of classical incom- 
pressible hydrodynamics for the vortex motion was first undertaken in [2]. The 
motion of a vortex half-ring in the presence of a half-sphere at an otherwise 
flat boundary was considered, which is solvable analytically because of maxi- 
mal symmetry by using a Legendre function of the second kind. The method 
is exact insofar as the (pinning) potential in which the vortex moves is exactly 
known in its relation to known geometrical quantities. This procedure involves 
no assumptions about the dynamics of the order parameter. The only funda- 
mental ingredient is the existence of the Magnus force acting on the vortex. 
This is essentially the approach which we will investigate in this paper. 

The approach of [3] consisted in evaluating transition probabilities between an 
initial many-particle wave function without and a flnal one with vortex. The 
many-particle wave functions were constructed from single particle solutions 
of the Gross-Pitaevskii equation [4], valid for a weakly interacting, dilute Bose 
gas. The main difference to [2], as discussed in [5], consisted in the existence 
of an (intermediate) deformed vortex ring state with a depression of the flow 
velocity out of the ring plane, leading to a smaller vortex energy. This also 
led to a different form of the tunnelling exponent by a logarithmic factor. 
Ultimately, though, this theory gave approximately the same value for the 
exponent as the collective co-ordinate approach of [2], provided one takes 
values of tunnelling parameters typically realized. 

In the recent past, there has been put some effort into numerical calculations 
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of instantons [6,7]. The instanton is a finite action solution to the Euchdean 
equations of motion and describes imaginary time motion under the tunnelhng 
barrier. In [7], an effective Lagrangian for a massive elastic string moving in 
a pinning potential near criticality was derived to calculate the shape of the 
instanton. The authors of [6], making use of the Gross-Pitaevskii equation, 
examined the nucleation of vortices in the streaming motion past an obsta- 
cle by solving the full nonlinear field equation in 1+ld. The intent there was 
to investigate the dynamics of the vortex linked to that of the condensate. 
With regard to this approach, it should be emphasised that use of the Gross- 
Pitaevskii equation for helium II on small (^) scales can at best claim phe- 
nomenological correctness, in this way comparable to density-phase functional 
theories [8]. 

Field-theoretical calculations, akin to those of pair creation in quantum elec- 
trodynamics, but for the corresponding one-dimensionally extended, closed 
object vortex ring in a U(l) superfiuid, were performed in [9-12]. The anti- 
symmetric gauge field tensor in the dual formulation of string dynamics [13] 
has been employed and the components of the 3-form field strength tensor 
identified with stringy generalizations of magnetic and electric fields. The the- 
ory is 'relativistic' in the sense that the propagation velocity of hght is replaced 
by that of sound in Lorentz invariant wave equations. The main problem of 
this formulation can be ascribed to the usage of an effective Nambu-Goldstone 
Lagrangian for sound on scales of O(^). This Lagrangian does not describe the 
real superffuid on these scales. In particular, because vortex velocities approach 
the speed of sound c^, these approaches give the hydrodynamic, i.e. large scale, 
mass of the vortex undue weight in the tunnelling exponent [14,15]. 

In summary, all these approaches have in common that they are mean-field 
theories, assuming that the quantum mechanical ffuctuations of the U(l) field 
4>{x,t) under consideration are on any scale relatively small compared to its 
expectation value =< (j){x, t) >, the order parameter. Such theories can not 
be derived in ^Ic from first principles on arbitrary length scales, because this 
dense superfiuid is a strongly coupled system, in which quantum fluctuations 
play a major role. On the atomic scale of the coherence length this implies, for 
example, that the quantum mechanical uncertainty for the density operator is 
of the same order as the local value of the density itself, the expectation value 
of this operator. Put very simple: One does not know (for purely quantum 
statistical reasons) if the particlc(s) in a coherence length sized volume is (are) 
there or not. The variable conjugate to density, the phase, has no definite value 
in such a small volume as well. In other words still, a collective central vortex 
co-ordinate of the quantum field ceases to be meaningful for the description 
of vortex motion on these small scales. 
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3 Semiclassical vortex quantum tunnelling as a geometric theory 

The aforementioned reservations about various (mainly field-theoretical) ap- 
proaches led the present author to a re- investigation [17] of the method [2] 

put forward already in 1972. The most appealing feature of this method is its 
simplicity and exactness in the semiclassical limit (cf. regions 1 and 11 in Fig. 
1). There is no reason for the description of the vortex as a stringlike entity in 
terms of a collective co-ordinate not to be a viable candidate for the investiga- 
tion of vortex quantum tunnelling, if the requirements for this formulation are 
carefully met. It should have become apparent from the above discussion that 
this method is in fact at present the only fully reliable one for the calculation 
of observable tunnelling probabilities in the actual, dense superfiuid helium 11. 

3. 1 The vortex action for a string 

To set up the general framework for the description of a vortex string in the 
hydrodynamic limit, it is convenient to use the dual formulation of string dy- 
namics in terms of the antisymmetric tensor gauge field h^^ already mentioned 
above in the field-theoretical context [16]. The appendix A contains a short 
summary of this formulation. The simplest way to get an intuitive picture of 
the dual formulation is that the tensor gauge field is a relativistic generaliza- 
tion of the stream function of classical hydrodynamics [18], i.e. it is chosen 
such that its (generalized) rotation is the velocity field of the fiuid. To begin 
with and introduce further into this formalism as far as we need it in what fol- 
lows, we then recall in addition that a semiclassical point vortex is analogous 
to a charged particle in a magnetic field. The stream function of classical hy- 
drodynamics plays the part of the negative scalar potential — a° = — $ acting 
on the vortex with a charge q = (mpoY^'^K. The vector potential Oj stems from 
the rotation of a 'magnetic' field which is antiparallel to the local circulation 
vector and related to the superfiuid density. These analoga can be readily gen- 
eralized to a stringlike entity by defining the above quantities locally on the 
string [16]. The Magnus force, a stringy generalization of the Lorentz force, is 
then derived in a Galilei invariant superfiuid from the variation of the action 

Sm = mpoK Jj^dtda (^^X'^ + kjX'^X') (1) 
= qj^dtda [-^ + aiX') 

with respect to the vortex co-ordinate X^{t,a) (we employ the convention of 
summation over equal indices). The line tangent is denoted X'* = dX^/da and 
the vortex velocity X^ = dX^/ dt. The bulk number density and the quantum 
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of circulation bear the symbols po and k, respectively. The parameter a labels 
points on the string. 

As expressed in the second line of equation (1), the first part of the action Sm 
is the generalization of the negative scalar potential in usual (non-relativistic) 
electro dynamic actions. The second part is analogous to the linear couphng 
term of the gauge vector potential = {mpoY^'^bijX'^ to the vortex velocity 
X*. The generalization from the point particle to the string essentially consists 
in introducing locality in a and admitting that the line tangent can point in 
any direction. The stream function becomes a vector with components ipi = h^i 
along the tangent vector. An isotropic gauge choice for the purely spatial part 
of h^y in Cartesian co-ordinates is given by hij — —{l/d)nijkX'', where d 
signifies the spatial dimension and riijk the unit antisymmetric symbol. This 
represents the isotropic solution of equation (A. 4). 

The local canonical momentum of a vortex per cr-length in a compressible 
superfluid 

= + = mpoK h,X'^ + V7 M^Xi (2) 



is a gauge dependent quantity. The quantity da is the measure of the 
string's proper length. The vortex momentum consists of two terms, of which 
the gauge dependent first one, Pf^^i is the contribution of the vortex- velocity 
dependent term in the Magnus action Sm (it corresponds to the electrodynam- 
ical minimal coupling term qai). The second one, the kinetic (vortex matter) 
part Pj"^™ = AfoATj, comes into play if the finite compressibility of the 
superfiuid is taken into account. This hydrodynamic contribution to the vor- 
tex momentum is acquired from the self interaction of the vortex string with 
phonons. The static mass of the vortex is in the present non-charged case 
Mo = Eseif/c^, i.e. its static self energy Egeif divided by the speed of sound 
squared [19]. In the given form the static mass is of relevance in the actual 
dense superfiuid if vortex velocities remain semiclassical, that is, much less 
than the Landau critical velocity related to the quantum many-body struc- 
ture of the superfluid. The hydrodynamic mass is, up to a factor of order unity, 
given by Mq ~ Mcore ln[P/(^e'^°)], where R is the infrared cut-off (distance to 
the next vortex or size of the vessel) and Co is a core constant, reflecting our 
ignorance of the core structure. The core mass Mcorc = '"^Po^r^^ is a measure 
of the mass per unit length of the normal fluid contained in the core. 

The contribution of P^^^ (coming mainly from region II in Fig. 1) is very small 
compared to that of the Magnus force acting in region I and II [14-17]. The 
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reason lies in the fact that the ratio of these two quantities, 



jkin 



1X1 



\X\c. 



(3) 



depends on \X\/cs as well as k,/{cs\X\) (= 0(^/|X|) in helium II). Both 
quantities arc necessarily <^ 1 if the vortex is to be described semiclassically. 
We used here that \bij\ is of order \X^\ in a Cartesian frame, and neglected 
the dependence of the self energy logarithm on the vortex co-ordinates. 

The vortex action has the canonical form 

T T 



Sv = J^dtdaX-P- JdtHv, 





(4) 



where the Hamiltonian Hy is the sum of static, kinetic, elastic and potential 
terms [16]: 



Hy — fda^/^ 



Q"] - mpoK I da (iv'c + ^u) ■ X' 



(5) 



Here, the vector Q is lying in the plane locally perpendicular to X' and mea- 
sures deviations from an equilibrium shape of the string. The stream function 
vector i/j (see Appendix A) is separated into its Coulomb part from the 
interaction of the vortex in question with other vortices and the part ij)^ stem- 
ming from the interaction with a (constant) background flow field. The first 
term in the square brackets represents the static, the second the kinetic, and 
the last one the elastic energy of the string. In this representation for i/y, we 
imposed the Coulomb gauge condition divP™'^ = 0, which results in the factor 
of 1/2 in front of the 'Coulomb' potential iI}q. 



3.2 Co-ordinates and momenta 



The rotational part of P™'^ satisfies rot P™^ — —mponX' (this follows from 
(A. 4)). Using a local righthanded co-ordinate basis e^, e^, = X' with 
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Fig. 2. The momentum f""^ defined from the Magnus force action (1) is a gauge 
dependent quantity. On every point a of the vortex Une it can point in an arbitrary 
direction of the local co-ordinate plane e^, eg. This direction is parameterized by 
the angle a{(j). 



determinant this yields the gauge invariant relation [17] 

d2Pr-diPr-rnp^i^Vg- (6) 



This equation expresses the well-known conjugateness of positions and mo- 
menta for a vortex line in an incompressible fluid (see, e.g., [20]). It appears 
here in its generalized form, valid for any choice of the co-ordinate basis on 
the string. Co-ordinates and momenta in different co-ordinate directions are 
no longer independent, just as for the analogous electrically charged particle in 
a very large magnetic field. The momentum direction depends on the local co- 
ordinate basis as well as on the particular gauge chosen. The Coulomb gauge 
for Hy in (5) is just one possibility. A convenient and often used gauge is to 
simply set = so that P^^"^ = mpon J dX^y^ (X). The co-ordinates and 
momenta can functionally depend on each other in a complicated, non-linear 
fashion, according to the choice of the co-ordinate system, which can be local 
in (7. Of importance in the context treated, though, is the fact that all the 
local momenta are given in terms of the local co-ordinates. 

A particularly well-known example [20] for conjugateness of the co-ordinates is 
provided by the Cartesian co-ordinates X,Y oi a. rectilinear line in z direction. 
In the gauge Py = 0, Px = mponY = hpoY and the commutator is 

[l,f]=.(27rpfV\ (7) 
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where pg is the 2d bulk number density in the x-y plane. For a circu- 
lar vortex in the gauge Pr = (cylindrical co-ordinates, y/g{X) = R), 
the canonical momentum per unit length of the azimuthal co-ordinate is 
P^(a) = {l/2)mpoKR^ = {l/2)hpoR^ and we have 

[Z,S]=ti27rp^^Y\ (8) 
where S is the operator of the surface area S = ttR"^ of the vortex ring. 

3.3 Volume elements in the tunnelling exponent 



The quantity of relevance for the description of tunnelling with constant energy 
E is the Legendre transform Se{E) of the Euclidean action S'e(T) (T is the 
period of motion in Euclidean time). The Euclidean action as a function of 
imaginary time is defined from the real time action by Se = —iS[t —it]. 



S,{E) = S,{T) -^T = S,{T) - ET 

= -i (j^da(j^ dX ■ P = (j^daj^ dK ■ P , (9) 

where we defined the imaginary co-ordinate K = —iX of the vortex[^ The 
meaning of Se{E) is most easily understood if we realize that the function 
exp[—Se{E)/h] is an energy dependent damping factor for a quantum mechan- 
ical wave function, traveling in imaginary time (with imaginary momentum), 
and penetrating a potential barrier (cf., in particular, the textbook [21] for a 
proper treatment of the semiclassical approximation). 

The peculiar property of a massless vortex, as seen in the last subsection, is 
provided by the fact that all the momenta have to be functions of the co- 
ordinates. If we forget about the small correction of the kinetic part to the 
momentum and set P = P™'^, using (6), we can thus bring the Euclideanized 
vortex action in units of Planck's quantum of action into the form[3] [17]: 

^ = Po /// V^dadZ'dZ' = pofiW (10) 



The co-ordinate differentials are represented by dZ^ = cos a dK^ + sinadK^ , 
dZ"^ = — sin a dX^ + cos a dX^ . The angle a{a) is the parameter giving the lo- 



K \s not to be confused with a wave vector. We could have chosen as well to 
incorporate the —i into the momentum P. Crucial is only that S^, is a real quantity. 
^ We omit the subscript V for the Euclidean action. 
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cal direction of the canonical momentum to be chosen subject to the constraint 
(6). The quantum of circulation was taken to he n = h/m. More generally, 
K — {Ny/Ns)h/m, with the vortex topological winding number and Ng the 
number of particles with mass m in the elementary boson. 

That the only variables in the problem are effectively co-ordinates yields an ex- 
pression for the tunnelling exponent which is in principle very simple. Given 
that we can solve for the motion of the local volume element on the line 
^du dZ^dZ^ along the whole path of the line in co-ordinate space as a func- 
tion of t, a, and its motion stays completely semiclassical, we can calculate 
the tunnelling exponent of the vortex to arbitrary precision. The motion of 
the volume element in (complex) configuration space (which represents simul- 
taneously the phase space) determines completely the semiclassical tunnelling 
exponent in the incompressible limit, that is, if this exponent is dominated by 
the contributions of region I in Fig. 1. This is always the case if -S'e(-E') » h. We 
will leave this semiclassical domain if the dimensionless action is approaching 
some number of order one or if characteristic variations of the co-ordinates 
take place on scales of the order of the core size ^. 

The integral (10) expresses the Bohr-Sommerfeld quantization of the number 
of particles A^*^''-* contained in the tunnelling volume O*^''-*: 

Se{E) = (A^('^) + j)h ^ Se{E)/h = 27r(iV('') +7). (11) 

The number 7 is of the order one and signifies the onset of the microscopic 
quantum regime. In the semiclassical approximation, N^'^'> ^ 7 must hold. 

3.4 Violation of Galilean Invariance 

At the absolute zero of temperature, a homogeneous non-relativistic super- 
fluid has Galilean invariance, that is, physical contents are invariant under 
co-ordinate transformations to any reference frame at constant velocity. If we 
approach absolute zero, which is what is actually realized in experiment, we 
expect the tunnelling rate to make no abrupt change as the temperature is 
lowered. Thus the result for the rate we obtain at T = should also be valid 
for temperatures slightly above T — 0. 

Because we can always transform to the rest frame comoving with the super- 
fluid, the tunnelling probability at T = equals zero if Galilean invariance 
remains unbroken: In the rest frame there is a tunnelling barrier of infinite 
height, the logarithmically diverging vortex self energy. Hence it is necessary 
to explicitly include the violation of Galilean invariance by a flow obstacle 
into any calculation of tunnelling rates for Galilean invariant superfluids at 
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absolute zero. 

The most likely location for the breakdown of invariance to happen in a pure 
superfiuid is the boundary of the superfluid having some surface roughness 
extending on scales much larger than ^. We will describe in the section which 
follows the simplest possibility to explicitly solve for the motion of a vortex 
in a boundary geometry of varying curvature. This will provide us with a 
picture of the geometrical restrictions on vortex quantum tunnelling in the 
semiclassical domain. 



4 Geometry Dependence of Vortex Quantum tunnelling 

It is obvious that to solve in general for the string motion of constant energy, 
according to the equations of motion resulting from (4) and the Hamiltonian 
(5), is a quite complicated task in some nontrivial boundary geometry. Deter- 
mining the tunnelling volume il^'^^ will be difficult in three spatial dimensions 
if we do not assume, like in [2], the highest possible symmetry of a unde- 
formed, massless ring vortex in the presence of a sphere. Even if we neglect 
any kinetic and elastic terms of vortex motion, the problem will require a 
quite formidable computational effort because to solve it, e.g., by the image 
technique will require a multitude of image vortices. 

Wc will thus cut the problem down on two dimensions to display the general 
properties and dependencies of the tunnelling volume on external geometry. 
This of course, then, does not describe any possible influence of waves trav- 
eling along the line. Their contribution to the Euclidean action, however, will 
be suppressed comparable to that of the kinetic term in (5), as long as semi- 
classicality is retained. 



4-1 Analytical hydrodynamic solution in two spatial dimensions 



The basic solution from which we start is that for a vortex in the presence of a 
half-circle of radius d at an otherwise fiat boundary (cf. Fig. 3). The complex 
plane of this original solution is called the Z-plane. The complex potential [18] 
generated by the image vortices and acting on the vortex at Zi is given by 




(12) 
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Fig. 3. The most simple nontrivial boundary problem solvable by the image tech- 
nique: A point vortex in the half space F > 0, which is filled with liquid, moving 
near a (half-) circle. The boundary conditions are satisfied by an image vortex at 
the plane and two image vortices of opposite strength inside the circle. 

There a three contributions, coming from the three image vortices depicted in 
Fig. 3. The first factor in the numerator stems from the image vortex at the 
plane F = with complex potential w^Z) = —i{K./2n) ln[Z — Zi], the second 
one is obtained by the circle theorem [18] as the image of the original vortex at 
the circle. Finally, the potential of the remaining -l-K-circulation vortex inside 
the circle, contributing in the denominator of the logarithm, completes the 
image vortex system, again by the circle theorem. 

The imaginary part of the complex potential gives the stream function ip — 
^[w], whereas the real part is the usual velocity potential. The first term in the 
denominator of the logarithm is incorporated into the static self energy of the 
vortex, itself = {mpoK^ /An) \n{/\Zi — Zi\/£^ , which is cut off by ^ and equal 
to half the energy of a vortex pair separated by \Zi^ Zi\. The expression for 
the potential in (5) is thus = -{h/2ti) In (1(^7^1 - Zi)/{(f/Zi - Zi)\\. 

In addition, we superimpose a constant external superflow of velocity u from 
right to left. This flow has complex potential 



w„(Z) = u[Z + 



(13) 



From (5), we can then infer the energy of the point vortex: 



E{Zi) = In 



[Z^ - Zi) [dyZ^ - Zi) 



^{dyz,-z^) 
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The energy is normalized by mpoii^/iTi, the characteristic energy of the 'par- 
ticle' vortex (corresponding in the electrostatic analogy to the charge squared 
divided by An [16]). We omitted in the above equation the contribution of vor- 
tex inertia, (l/2)Mo|Zip (where Mq — Eseif/c^), i.e. wrote down the energy 
in the incompressible limit. 

We wish to map by a conformal transformation the region outside a boundary 
surface with varying curvature radius, lying in the 2;-plane, to the domain out- 
side the circle. Any such transformation can be written z = aoZ + J2'i^=o bnZ~"', 
where ao, bn are some coefficients and Z — dexp(ix) is on the circle [18]. We 
would like to invert this relation to obtain the solution for the boundary surface 
directly from that for the circle. The simplest possibility to do so is supplied 
by choosing 60 = 0, 6„ = for n > 1 and scaling ao to unity, which leads to 
the celebrated Joukowski Transformation 

z^Z-f/AZ, (15) 



which maps the outside of the ellipse with half axes a,b {a < b) to the outside 
of the circle of radius d = {a + b)/2. The parameter / is defined by P — 
b'^ — a'^. The inversion of the transformation (15) 2Z = z-\-\/z^ + P gives, upon 
insertion into the solution for the circle, the expression for the vortex energy 
in the presence of the ellipse. It can be seen from the structure of the inverse 
transformation that the use of elliptic co-ordinates, defined hy z — I sinh is 
most convenient. Here, = % + iry is the new complex co-ordinate, related to 
a Cartesian system by a; = / sinh % cos r/, y = / cosh % sin r/. The co-ordinate 
lines are confocal ellipses and hyperbolas, and 2Z = /exp(C). 

In what follows, we consider the solution (14) in the half-plane y > 0. The 
normalized energy as a function of the elliptic vortex co-ordinates Xi, 771 takes 
the form 



b exp(xi - Xo)| sin?7i| sinh(xi - Xo) 



(sinh^(Xi - Xo) + sin^ 771)1/2 
Xo) I sin 771 1 , 



inula + b) , . 

- smh(xi 



(16) 



where Xo = artanh(a/6) is the co-ordinate specifying the ellipse surface. The 
shape of the energy in the X'V plane is shown in Fig. 4. 



The velocity u will henceforth be scaled with the characteristic velocity vl — 
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Fig. 4. Shape of the potential barrier (16) with the choice xo = 0.175. The corre- 
sponding elhpse with 6/a ~ 5.7 is shown in Fig. 5. The (large) velocity u = 0.08 is 
in units of vj^ = k/2ttE^ and the ratios a/^ = 2, s/^ = 1. The zero of this normalised 
potential energy is shifted by Eq, defined in (17). 

«;/27r^, expressing a velocity limit for the validity of the semiclassical approx- 
imation. (The Landau critical velocity of roton creation ~59m/satp~l 
bar equals if ^ ~ 2.7 A.) The velocity u — 0.08 chosen for the potential 
in Fig. 4 is large and near critical in the sense that the local velocity at the 
ellipse top is near v^. It appears from this graph in real co-ordinate space that 
the vortex will predominantly tunnel in the x-direction around the ellipse top 
(?] ^ 7r/2), because there the potential barrier is shallowest. 



4.2 Geometric restrictions for semiclassicality 

Of course it is possible to invent more intricate conformal transformations 
to map the circle solution onto fancy shapes in the 2;-plane. But the ellipse 
already possesses the feature crucial for an analysis of vortex motion in a 
geometry deviating from the highest symmetry of Fig. 3: It has curvature 
varying continuously between a/b'^ in the x- and b/a^ in the y-direction. This 
leads to a restriction on the possible vortex paths near the boundary [17], 
visualized in Fig. 5. 
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Fig. 5. Two vortex paths of constant energy near the ehipse. Whereas the vortex on 
path 1 with approximately zero energy, E 0, does not manage to pass by without 
coming closer than ^, the second one, having energy E = Eq, defined in (17), is able 
to do so. The velocity u is to be understood 'at infinity'. 



On the first path we consider a vortex having approximately zero energy 
^ 0, starting far away from the ellipse with a distance of the order ^ 
to the boundary. As it approaches the ellipse and tries to get around it, it will 
inevitably hit the boundary because it does not have enough energy to cross 
the sharp ellipse top at a safe distance. 'Safe', that is, such that the semiclas- 
sical approximation can be applied. For that purpose, the vortex has to have a 
distance to the boundary exceeding the many-body core size ^. This imprints 
a restriction on the value of the possible tunnelling energies for a geometry in 
which the curvature radius decreases (the curvature increases) while the vortex 
is passing by: The energy has to be chosen with a value at least high enough 
such that the vortex is able to pass the complete boundary at a distance larger 
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than ^. Introducing the closest approach distance s to the boundary (see Fig. 
5), and observing that for Sx = Xi ~ Xo ^ the distance of the vortex to the 
elhpse top is given by 6b ~ a6x, we impose Sxivi = vr/2) = s/a, where we 
assume b ^ a ^ s and s > 0(^).[3 Under these conditions, the vortex center 
on path 1 in Fig. 5 will pass the top at the distance 6b ~ {a/b)^, whereas the 
path 2 vortex with energy Eq = Eo{a, b, s, ^, u), 



En 



~^ - exp 6x tanh 6x 




^ smh 6x (17) 



vl \ a / ? ' 

passes at 56 ^ s. The energy Eq is the energy needed by the vortex to remain 
completely describable in semiclassical terms on its way along the ellipse. In 
general, it depends on u, but in the low velocity limit 2u/vi <^ 1, it is given 
by Eq ~ ln[(6/a)(s/^)]. The effect of this nonzero energy can be interpreted 
as the rescaling of ^ — > ^ exp £"0 in (16). The energy is zero with respect to this 
rescaled ^. For small m, the core size rescales as ^ — {b/a)s. 

That the validity of the semiclassical approach enforces that we introduce 
another geometrical quantity, s, is a restriction of quite general character. 
It is of relevance for any attempt to describe tunnelling semiclassically in a 
realistic, non-spherical geometry, i.e. when the boundary and thus the path of 
the tunnelling object near it is not of S"' symmetry. A semiclassical description 
is valid only if the quantum core structure of the tunnelling object is not 
touched upon during its motion along the boundary. A pinning potential for 
the vortex moving in the superfiuid stems in general from some flow obstacle, 
in our case the ellipse. Any phenomenological ansatz for a pinning potential 
usually employed in tunnelling calculations which has curvature perpendicular 
to the applied flow larger than parallel to the flow will have to take into account 
that the object can approach the surface within its core size, invalidating the 
collective co-ordinate description. 

It is to be noted that the energy is in units of mpoK'^/Air ^ 0.82 K/A (at p ^ 1 
bar and in three dimensions). For realistic values of the parameters in (17), 
the values of Eq cover the same range as the phonon-maxon-roton spectrum. 
Energywise, the trapped vortex thus can not be distinguished from an elemen- 
tary excitation. It could have been excited thermally and remained trapped 



^ In the relation s > 0(.^), the quantity O(^) means 'a value in the order of ^ by 
definition', as there is no sharp distinction between the 'inside' and 'outside' of a 
quantum vortex. Additionally, the choice for the lower limit value of the distance s in 
units of ^ depends on the value of the core constant Co = 0(1), parameterizing the 
many-body core structure in the vortex energy logarithm ln[R/ (^ exp Co)] (exp Cq = 
1 for the point vortex considered, which has R = 2Y). 
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in a pinning center during the cool-down of the superfluid to very low temper- 
atures. From this and the above analysis it follows that it is scmantically not 
appropriate to use the term 'nucleation' if we remain completely semiclassical 
and define it to mean 'creation from nothing pre-existent', that is, from the 
zero of energy. Experimentally, it will be impossible to distinguish the tun- 
nelling of a pre-existing small energy vortex at a rough boundary from the 
true nucleation event of a nascent vortex there, if no direct means to control 
the microscopic dynamics can be provided. 



4-3 The tunnelling area 

In the effectively one- dimensional problem we are considering [X is the loca- 
tion of the vortex, Px — mp^KY — hpoY its momentum), there always exists 
a closed vortex path in phase space. For the tunnelling trajectory this path 
is in the K-Y plane and the closed vortex path is running underneath the 
barrier of Fig. 4. It has in the limit of small velocities the shape shown in Fig. 
6. 

V 




K=-iX K^=^Y, 

Fig. 6. The closed vortex path giving the action (18) in the low velocity limit. The 
first part Ye corresponds to the analytically continued path 2 of Fig. 5 along the 
ellipse surface. The second part of the closed path, Yn ~ Yq, obtained from the 
second solution of (16) for the constant energy in (17), represents the border line 
to a free vortex. 
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The action (9) in the form of (10), which follows from the closed vortex path 
in the low velocity limit [17], 



Krr 







a s\ j K, 
\ b a) \ Anu 



In 



K, a/b 
2t:u s 



(18) 



is bounded from below by the fact that the tunnelling can not be 

arbitrarily reduced without violating semiclassicality, i.e. by the introduction 
of s. Geometrically, the curvature radius of the vortex path near the flow 
obstacle can not be made arbitrarily small to fit a given shape of the asperity 
on the boundary, so that the eccentricity e ~ 1 — 0.5/3^ of the path can not 
approach unity arbitrarily close. The bigger of the two quantities a/b,s/a > 
^/a decides about the lower limit of ^2^^^ attainable in semiclassical tunnelling. 



The materialization distance Yq of the vortex is in the low velocity limit given 
by (in units convenient for the discussion of experiments which follows below) : 



Y. 



0.8 nm 



' vl a/b^ 

w[10m/sec] \u s/^^ 



In 



(19) 



In three dimensions, the radius of the materialized vortex half-ring, Rq — 
{k/2ttu) ln(- • •) = (1.6nm/u[10m/sec]) ln(- • •), bears an additional factor of 2 
in front of the logarithm. 



5 Prefactor estimations; Relation to experimental data 



The tunnelling probability has the form P{E) = Aexp[—Se{E)/h], where A is 
the so-called prefactor depending on the fluctuations of the quantum variables 
around their classical value. Apart from the considerable difficulties in evaluat- 
ing prefactors in general, an accurate calculation of A in a dense superfluid like 
helium 11 is in principle not possible due to the lack of a microscopic theory. 
It is, however, feasible to get an idea about the value of this prefactor within 
about two orders of magnitude. This is all what we need, because, as we will 
see below, the variations of the tunneUing exponent Se{E) » h with the (geo- 
metrical) parameters of the problem will dominate anything which is actually 
observable (provided of course that a semiclassical, large scale description is 
appropriate). The simplest possible idea about the prefactor is gained by con- 
sidering the frequency cUa of a particle oscillating in a metastable well. Then, 
neglecting the influence of dissipation on the vortex motion, within about one 
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order of magnitude A ~ [27]. In the thermal activation hmit, i.e. in the 
Arrhenius law case P = (c<Ja/27r) cxp[— [//A-^T], Ua is the frequency of oscil- 
lations in the metastable well against a barrier of height U. The frequency 
= '^a/27r can thus be generally understood as a measure of the number of 
times per second the vortex bounces against the potential barrier, trying to 
get free. Wc have no possibility to describe the vortex state (at the boundary) 
quantitatively, but we are able to conclude on the order of magnitude of Ua 
by taking into account that there exists a surface layer of vorticity of width 
^. Because the superfluid density goes to zero at the boundary and heals back 
within ^, the energy needed for the activation of vortices vanishes within this 
distance [3] . The frequency of motion of these vortices should then be of order 



a;o = ^ = 4.87 • 10^ ( ^ ) , (20) 




which is the cyclotron frequency of vortex motion. We scaled ^ with the 
Lennard- Jones parameter ulj = 2.556 A of the ^ie atomic interaction. Ap- 
proximately the same estimate is obtained if we directly take into account 
compressibility in the form of a finite speed of sound. The rigidity (spring 
constant) of the vortex against deformations will scale as 1/^^, its mass as 
1/c^, so that 

^, = ^~9.23-10iW^^, (21) 



where we took the pressure to be p ~ 1 bar. It should be stressed that the 
approximate coincidence of the estimates (20) and (21) is a particular feature 
of helium II, where k ~ 2^0^. If we adhere to the description of a massive 
vortex of size ^ moving in some (regularized) potential, the equation (21) is 
to be used, whereas (20) is the natural vortex frequency associated with the 
scale ^ alone. 

Assuming that the prefactor can vary between A ^ 10^° • ■ ■ 10^^ Hz, its loga- 
rithm In A ^ 23 ■ ■ -28. It is thus obvious that variations of the action in the 
tunnelling exponent with geometrical factors will outweigh by far any varia- 
tions in the prefactor. The tunnelling exponent in (11) involves N, the number 
of particles in the tunnelling volume. Corresponding to the conceivable values 
of the prefactor, it will have to vary between iV = 4 • • • 6, say, for tunnelling 
events to be observable within a reasonable span of experimental time. Though 
this number is small, it is not impossibly so for the hydrodynamic treatment to 
make sense at least asymptotically {asymptotic to the microscopic nucleation 
regime) . 
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The experiments in which one hopes to observe intrinsic^ vortex quantum 
tunneUing at irregular boundaries [22] are using flow orifices of sub-micron size 
[23] . An oscillating superflow of frequency in the order of Hz is driven through 
the orifice. It is observed that at a certain critical value of the amplitude of the 
diaphragm driving the flow through the hole, there is an instantaneous (on the 
scale of the driving frequency) breakdown of the amplitude which is quantized 
[24] . This quantized dissipation event is associated with a vortex generated at 
the orifice walls, subsequently crossing all the streamlines of the flow, thereby 
causing a phase slip event [25], which draws a quantized amount of energy 
from the flow. The critical velocity, at which the phase slip takes place, is first 
linearly increasing with decreasing temperature and then saturates at T ~ 150 
mK to a plateau (/)Iat within experimental resolution). This plateau has been 
conventionally associated with the quantum tunnelling of half-ring vortices 
through a potential barrier [26]. 

In the experiments performed so far, the mean critical velocity through the 
orifice on the plateau is of the order 5-10 m/s. Though this is not necessarily 
equal to the velocity u at 'infinity' (think, e.g., of a sharp spike on a smoother 
asperity at the orifice wall), the radius Rq of the materialized vortex half- ring 
can be estimated to be in the order of nanometers. This should be sufficiently 
big for the hydrodynamic approach to make sense. The crossover temperature 
To from thermal activation to quantum tunnelling is, for negligible dissipation 
and small u, given by 27ikBTo = hub, where Ub is the oscillation frequency 
of vortex motion in the inverted potential (see, e.g., [27], [28]), which yields 
uja — uJb = 8.2 ■ IO^^Tq [100 mK] Hz. This is about an order of magnitude 
smaller than the estimates (20) respectively (21). There are, however, several 
facts making a direct comparison of these values questionable. First of all, these 
estimates can give only a rough idea about the true dynamical behaviour of 
a many-body vortex near the boundary. It is conceivable, for example, that 
the effective 'spring constant' of the vortex against deformations is lowered 
compared to the semiclassical estimate in (21) because of the many-body 
quantum uncertainty of its position. In addition, we have seen in the preceding 
section that a vortex will in general not be in a state with zero energy at the 
wall. It will rather be in a state with the energy Eq in (17). This can also cause 
a change in the prefactor. Furthermore, the prefactor is in general a function 
of driving velocity u and temperature T [28] and at the measured crossover 
temperature in terms of the critical velocity not necessarily equal to its value 
at zero temperature. This value will presumably be closer to the estimates 
(20) and (21). 



^ The term intrinsic is meant to be the opposite of 'supported and supplied by 
pre-existing vortices of macroscopic size', i.e. extrinsic. 
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6 Conclusion 



There is no theory of vortex quantum nucleation in the dense Bose super- 
fluid hehum II which could resolve the problems connected with theoretical 
issues and experimental facts. There is, though, a consistent theory of vortex 
quantum tunnelling in the semiclassical domain, which we represented here in 
its formal requirements and geometric implications. It is again to be empha- 
sised that this theory does not make any assumptions about the dynamical 
behaviour of the superfluid on scales approaching the microscopic one. This 
implies on one hand that it necessarily can not supply a theory of vortex nu- 
cleation. On the other hand, it gives semiclassical bounds on tunnelling rates 
which can be experimentally observed. It should have become apparent from 
our discussion that to measure the consequences of the actual microscopic nu- 
cleation event and, in particular, influence the (microscopic) parameters which 
govern its probability, will be a task quite difficult to realize. If the tunnelling 
event takes place on sufficiently large semiclassical length and small velocity 
scales, i.e. is sufficiently semiclassical, the only remnant of microscopic dy- 
namics we can expect to be involved in the tunnelling rate is the core size ^. 
But this is only appearing through a logarithm, which means that its influ- 
ence on the tunnelling rate is quite subordinate as compared to the tunnelling 
exponent variation with the path geometry of the vortex, imprinted on it by 
the flow obstacle. 

The very fact of quantum tunnelling at boundaries needs further proof so that 
the predictions of tunnelling theory can be compared to that of classical insta- 
bility mechanisms (investigations in this direction are found under reference 
[29]). One such proof could consist in the comparison of critical velocities for 
chemically identical oriflces of equal global sizes, which have different micro- 
scopic surface structures. If the result of such measurements is negative, i.e. 
there is no reproducible difference in critical velocities, there is no quantum 
process taking place describable by (semiclassical, hydrodynamic) means of 
tunnelling at irregular boundaries. The experimental information at present 
is too sparse to give a flnal and conclusive answer. 
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Appendix 

A The dual formulation 

In this appendix, we introduce the reader to the basics of the dual formu- 
lation of string dynamics. For ease of reading the formulas are given in the 
conventional index notation, which will be familiar to the majority of the 
readership, as well as in a co-ordinate independent formulation using p-forms, 
which shows most clearly their geometrical meaning [30]. We give the covari- 
ant formulae first, subsequently reducing them to their Galilean limits, which 
are of relevance in superfiuids like helium II. 

In the hydrodynamic limit, a neutral superfluid of spinless bosons is described 
by an order parameter function = p^^'^{x) ex.Y>\i6{x)] of the spacetime point 
X = {t,x) (the velocity of light c = 1). Its absolute square is the density p{x) 
and its phase 9{x) the velocity potential of the fluid. The fundamental conser- 
vation law derived from this semiclassical U(l)-description of the superfluid is 
that of the hydrodynamic mass four-current J = mpv = hpdO = hpd^O dx^ ^ 
where 9^ = djdx^ . The most familiar mathematical form of mass conservation 
is provided by dtP + div(pw) = 0. This is expressed covariantly as 

V^J'' = ^ dA*J = (<Z=cZx"9„). (A.l) 

That is, the dual of the current J = J^dx^ is a 3- form *J — J'^e^^.ajs dx'^Adx°'A 
dx^, which is closed. This statement is equivalent to a vanishing covariant 
derivative V^J^ of the current vector with components J^. We define the 
field .strength H by *J = nipoH, where po is the bulk constant density. The 
field H^i^a is totally antisymmetric in its three indices and has, by definition, 
only four independent components. In a simply connected region, H is exact, 
i.e. it is the exterior derivative of a gauge 2-form b = h^^dx^^ A dx'^: The 
field strength H is invariant under gauge transformations h b + d A A, 
where A = A^dx" is an arbitrary 1-form. In components: The replacement 
bj^v b^y + d^Ai, — dyA.^ leaves H^^a unchanged. We have 

*J = mp*v = hp *d9 = mpo d Ab = mpoH, (A. 2) 

Hapj = d^bp^ + d^bafs + dph^a = {p/ Po)v''e^ai3^ , (A.3) 

which defines the relation of p, 9 and b. The dual transformation in the con- 
ventional sense [13] is obtained if p — po, i.e. if we consider compressibility 
to be negligible in the domain of interest. Then, b has only one degree of 
freedom corresponding to the order parameter phase 0. For our purposes, this 
reduction to a single degree of freedom amounts to not considering the cores 
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of vortices. This is reasonable, because in the dense superfluid hehum II the 
semiclassical U(l) description is not vahd in the vortex core. 

The superfluid we are deahng with is only Galilei invariant, not Lorentz in- 
variant, so that by writing (A. 3) with p = po for spatial and temporal indices 
separately and taking the Galilean limit, we get the relations 

- \^nijk = dkbij + djbki + dibjk = Hijk , (A.4) 
y/guijk v' = dkipj - djTpk , (A.5) 

where ipi = boi is the vectorial version of the stream function of classical 
hydrodynamics [18]. The determinant of the spatial co-ordinate system we are 
using is designated g and riijk — n[ijk] — ±1 is the unit antisymmetric symbol. 
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